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Global approximants with renormalization scale invariance in pQCD 
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Truncated perturbative series (TPS’s) of any observable have the unphysical dependence on the choice of the 
renormalization scale (RScl). The diagonal Pade approximants (dPA’s) to any TPS of an observable possess the 
favorable property of being invariant in the large-/?o limit. This means that they are invariant under the change 
of the RScl fj, 2 when the “running” coupling parameter a a (g 2 ) evolves according to the one-loop renormalization 
group equation. We present a method which generalizes this result - the resulting new approximants are fully 
RScl-invariant in the perturbative QCD (pQCD). Further, we present some numerical examples. Both the dPA’s 
and the new approximants are global, i.e., their structure goes beyond the usual (polynomial) TPS form and thus 
they could reveal some non-perturbative effects. 


The contribution is based partly on [Q. It 
contains additionally some numerical examples. 
Since the construction of our approximants is 
related with the diagonal Pade approximants 
(dPA’s), we first explain what PA’s are. 

1. What is Pade approximant ? 

Suppose we have a physical quantity F{z) 
which depends on the parameter z. In the ex¬ 
pansion of F(z ) in powers of 0 

F(Z) = f 0 + flZ + ---f n Z n + ---, (1) 

the fj' s are in principle calculable. Suppose 
that /o,...,/ l~\-m have been calculated, i.e., the 
truncated perturbation series (TPS) Fy L+M ^ (z ) is 
known. Then, the PA [. L/M]p(z ) of order L/M to 
F(z) is defined as the ratio of two polynomials, 
of degree L (nominator) and M (denominator), 
such that, when expanded back in powers of z, it 
reproduces the first L+M+l terms of (|Ij). In gen¬ 
eral, this condition determines uniquely the PA. 
It is the minimal condition that any approximant 
to the TPS F[ L+M ](z) has to fulfill. PA has, in 
addition, the favorable property that it goes be¬ 
yond the analytic form, showing structures (pole 
singularities in the complex z plane) not explic¬ 
itly contained in series (1). We call such approxi¬ 
mants global , since they can give us clues to some 
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nonperturbative properties of F(z) at large |^|. 
The diagonal PA’s (dPA’s) are those with L = M. 

2. Construction of new approximants 

A generic observable S in pQCD is 

S = a(q 2 )f(q 2 ) = a(q 2 )[ 1+ 

r 1 (q 2 )a(q 2 ) + --- + r n (q 2 )a n (q 2 ) + ■■■], (2) 

where q 2 is a chosen renormalization scale (RScl), 
a(q 2 ) = a s (q 2 )/n. S is RScl-independent. How¬ 
ever, available is only a TPS S\ n i 

S[ n ]{q 2 ) = a(q 2 )f [n] (q 2 ) = a(q 2 )[l+ 

r 1 (q 2 )a(q 2 ) + --- + r n (q 2 )a n (q 2 )}. (3) 

It has unphysical RScl-dependence (truncation). 

How to extract as much information as possible 
from the available TPS (||)? Firstly, any approxi¬ 
mant should fulfill the minimal condition (cf. pre¬ 
vious Sec.). Secondly, the full observable S con¬ 
tains, in general, non-perturbative effects not ex¬ 
plicitly manifested in power series ([|) - this leads 
us to consider global approximants (cf. previous 
Sec.). Thirdly, S is RScl-independent, so it is nat¬ 
ural to expect that RScl-invariant approximants 
bring us closer to S. So, the question here is: How 
to construct global approximants which are based 
on the TPS S[ n ](q 2 ) and are RScl-independent? 

It turns out that a partial answer to this ques¬ 
tion is the diagonal Pade approximant (dPA) 
[M/M\s(a) @. To see this, we recall that the 



2 


evolution of a(p 2 ) with the change of RScl p 2 is 
governed in pQCD by the RGE 

~ 7 ~j d ( a 2 , = ~Poa 2 {l + cia + c 2 a 2 H-). (4) 

dm{—p z ) 


In the large-/3o limit (i.e., ci = C 2 = • • • = 0) 
a { 1 L) (p 2 ) = a(q 2 ) [l + \n(p 2 /q 2 )/3 0 a(q 2 )\ 1 . (5) 

Thus, the RScl change q 2 1 —> p 2 results in the 
change of z = a(q 2 ): z i—> z/(l + bz). The 
dPA \M/M]s{z) is invariant under this argu¬ 
ment transformation, i.e., RScl-invariant when 
a(p 2 ) evolves according to the one-loop RGE. 
[M/M]s(z) reproduces the TPS S'pM-i] © to all 
the available powers (cf. also the previous Sec.). 

Is it possible to go beyond this large-/?o ap¬ 
proximation? The question was raised first by 
the authors of © The answer: it is. 

The dPA approach has to be extended. We 
describe the algorithm leading to the new, fully 
RScl-invariant, approximants to S. Introduce 


k m (a(q 2 )) 
km(a) 


1 d m a(p 2 ) 


a(q 2 ) m! d(lnp 2 ) m p 2 -q 2 ’ 
(—l) m /3™a m [l + 0 {cia)\ , 


( 6 ) 

( 7 ) 


where higher orders in © can be calculated by 
RGE 0). We rearrange power series (© for S 
into perturbation series in k m (~a m ) 


S = a(q 2 )f(q 2 ) = a(q 2 )[ 1 + fi(q 2 )k 1 (a(q 2 ))+ 

H-1 ~ fn{q 2 )k n (a(q 2 )) + ■ ■ ■], (8) 


where f n is uniquely determined by the first n 
coefficients r, of (0). For example, fi=—ri/Po, 
/ 2 =(—Ciri+r2)//3o) etc. Thus, knowing the TPS 
S[2M-1] (q 2 ) (n,- ■ ■ ,r 2 M-i) in powers of a, we 
know the corresponding TPS in k m ’s of © up 
to (and including) the JW-i-term. Define the 
(power) TPS, by the large-/?o substitution [cf. (0)1 

a(q 2 )T [2M ~ 1] (q 2 ) = o[l - /i/? 0 a+ 

+ • • • + f 2 M-i(-l) 2 M - 1 p 2 0 M ~ 1 a 2M - 1 }. (9) 

The dPA \M/M](a) for this (power) TPS can be 
decomposed into a sum of simple fractions 

M / 2 \ 

= 1°) 


We now denote 

p 2 = q 2 exp [Uj(q 2 )}, i.e. Uj(q 2 ) = In (p 2 /q 2 ), (11) 
and therefore, in view of ©, rewrite © as 
M 

[M/MU = J 2 ^ 1 L \p 2 ) . (12) 

3 =1 

Replace here a^ ll '\p 2 ) 1 —> a(p 2 ), where a(p 2 ) is 
evolved from a(q 2 ) via the full RGE © 

M 

4 M/M] = • ( 13 ) 

3 = 1 

This approximant satisfies all the conditions that 
we set forth at the outset: a) when expanded back 
in powers of a(q 2 ) (q 2 is the original RScl), it re¬ 
produces all the terms of the TPS 5[ 2 m-i] (<? 2 ) 
(p|); b) it is global since it is a modification of the 
dPA approach (the latter is global); c) it is fully 
RScl-invariant in the pQCD-sense, i.e., indepen¬ 
dent of the initial choice of the RScl q 2 , where 
a{q 2 ) evolves according to the full available per¬ 
turbative RGE ©). Explicit proof is given in JT|] 
(first entry). In fact, coefficients ciq and squared 
momenta p 2 are all RScl-invariant. 

Approximants ([© can be applied directly only 
to TPS’s S[ n ] © with an odd number n = 2M— 1. 
In pQCD, jS[;l] are available for many observables, 
S'[ 2 ] for a few, S pj for none. Thus, (111) is appli¬ 
cable only to TPS’s Spj. In this case (M = 1), 
( 111 ) reduces to the effective charge (cf. [Q). To 
apply 0) to Spi^ 2 ), a modification is needed. 
Introduce a new observable S = S * S 

S = ( S ) 2 = a(q 2 )F(q 2 ) = a(q 2 )[ 0 + la(g 2 )+ 

+R 2 {q 2 )a 2 (q 2 ) + R 3 (q 2 )a 3 (q 2 ) -|-], (14) 

where l?i = 1, R 2 = 2rq, R 3 = r\ + 2r 2 , etc. Hav¬ 
ing S[ 2 ](q 2 ) (|), we have %](g 2 ) (i.e., up to and 
including the l? 3 -term). We apply 0) to S 3 
( M = 2). Although the leading order term in ( |l4| ) 
is 0 * a(q 2 ), the algorithm survives, results in 

S = Af 2] + 0(a 5 ), S = yjA [2/2] + 0(a 4 ). (15) 
\J a § /2] =yja 1 HpI) - a{pl)\ , (16) 

where aq and p 2 are determined by rq, r 2 and 
/?o,Ci,c 2 (cf. § , last entry). Approximant © 
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is fully RScl-invariant. The scales p 2 may be com¬ 
plex in some cases, but the result is real. 


3. Specific numerical examples 


3.1. A case of a Euclidean observable 

Consider the Bjorken polarized sum rule 


dx 


g[ P \x,Q' 2 ) - g[ n) (x,Q' 2 ) 


( n )i 


[\g A \/(3\g v \)][l-S(-Q 2 )], 


(17) 


p 2 = —Q 2 < 0 is 7 * momentum transfer. TPS 
S[ 2 ](—Q 2 ) is now known: ry = 3.584 (j^|); V 2 = 
20.212 (§). These ry’s are in the MS scheme, 
rif = 3, and the chosen RScl g^ Scl = —Q 2 - We 
take 1.76 GeV [a{-Q 2 ) =0.083], 

First we note that any approximant can be used 
also to predict the next coefficient r 3 , by reex¬ 
panding the approximant back in powers of a. We 
give in Table 1 results of various approximants to 
5 of ©: PA’s, approximant (pb|), and the effec¬ 
tive charge method (ECH Q, we set C 3 ECH ^ =0). 
In brackets, results are given for another scheme 
(c 2 = 3, C 3 = 0). We used the known C 3 parameter 
of the MS scheme [|7j wherever possible, i.e., MS 
is characterized by C 2 =4.471 and 03 = 20.99. The 
results of method © differ somewhat from those 
of other methods. Further, values of S(—Q 2 ) 
(cf. second column) predicted by global methods 
show up some scheme (C 2 — )dependence (~ 1 %). 
The RScl-dependence is another source of ambi¬ 
guity in the results of the PA’s. For example, if 
changing RScl Q 2Q, the result of the [1/2] 
PA changes by 3.5%, that of [2/2 ] 1 / 2 by 2%, and 
the original TPS Spj by 11%. There is no Rscl- 
dependence for approximants ([L 6 ]). The result 
S(—Q 2 ) of the ECH method is scheme- and RScl- 
independent, since this method is local (i.e., a 
specific choice of RScl and scheme). 

If changing the scheme more drastically, e.g. 
MS 1 —> ‘t Hooft scheme (02 = 4.471 1 —>• 0; 03 = 
20.99 1 —> 0), predictions for S(—Q 2 ) of the global 
methods [ 1 / 2 ], [ 2 / 2] 1 / 2 and ( |l 6 |) change by 1 . 2 % 
(0.1273 0.1260), 1.2% (0.1326 0.1342) and 

4.4% (0.1378 1 —> 0.1.439), respectively. Scheme 
dependence (i.e., C 2 -dependence) of © does 
present a problem in this case. 


3.2. A case of a Minkowskian observable 


R(p 2 = s) = 


cr(e H 


7 


hadr.) 


cr(e _t 


p+p ) 


(18) 


Here, p 2 =(p e +pg) 2 =s>0 is Minkowskian. Ap¬ 
proximants generally lose predictability when ap¬ 
plied to such observables ||. The associated 
Adler function D(p 2 =—Q 2 ) (Q 2 > 0), RScl- and 
scheme-invariant, is Euclidean 
/»00 

D(-Q 2 )=Q 2 dsR(s)(s + Q 2 )~ 2 . (19) 

Jo 

We apply approximants to it. The TPS up to 
NNLO has been calculated for D [j|. For n/=5, 
^RSci = Q 2 ■> in MS, it is 

D = y[l +d {iul \ d (full > =d + d { - u \ (20) 

d(-Q 2 ) = a(-Q 2 )[l + 1.4092a(-Q 2 ) 

-0.6812a 2 (—Q 2 ) + •••], ( 21 ) 

d (u.)(_Q 2 ) = a 3 (—Q 2 ) (-0.3756)+ •••. (22) 


We take yjcf 2 = 34 GeV, a(-Q 2 ) = 0.0452. The 
d (u.) 

is from light-to-light diagrams, should not 
be included in the approximants since the re- 
sunnnations cannot “see” separately diagrams of 
fundamentally different topologies. Reexpanding 
the approximants to d(—Q 2 ) in powers of 
a(— Q 2 ), we predict ds in series (| 2 l|) (note: di = 
1.4092, d 2 = —0.6812). Then relation © allows 
us to obtain the coefficients 77 of the expansion 
of r(s) in powers of a(—s) [R= (ll/3)(l+r( full i), 
r(fuU) = r -\-r( 1 - 1 -), r ( l L ) =d^ LL ), (— Q 2 = — s)] 

r(s) = o(-s)[l-t-rifl(-s)-t-], (23) 

with 77 =di = 1.4092, r 2 =d 2 — Tr 2 pQ/3= —12.767, 
r 3 =^ 3 —89.190, etc. Thus, we can predict also 
the coefficient r 3 of r(s). Results are given in 
Table 2. In brackets are results in the ’t Hooft 
scheme (Cj = 0 for j > 2). Predictions of the 
PA methods and the ECH are clustered. Those 
of our method are slightly, but significantly, de¬ 
tached from them. Scheme (c 2 -)dependence for 
the global approximants is weak (cf. last two dig¬ 
its in 2 nd and 4 th columns). Further, if the RScl 
is changed 9 rs c i = — Q 2 ^ Q 2 /4, the d(—Q 2 )’s of 
the PA’s change about twice as strongly as when 
the scheme is changed C 2 = 1.475 1 —> 0. Approxi¬ 
mant © does not change when RScl changes. 
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Table 1 


Predictions of various approximants for the Bjorken polarized sum rule (with nf = 3) in the MS (c 2 = 4.471) 
and the c 2 = 3 scheme (in brackets). \J~Q 2 = 1.76 GeV; RScl chosen: g^ gcl = —Q 2 ; a(—Q 2 ) = 0.083. 


Approx. 

S(-Q z ) pT - 

'3 


% (TPS) 

0.1192 (0.1175) 


0.1192 (0.1175) 

ms 

v 7 [2/2] s 2 

W' 

0.1273 (0.1261) 
0.1326 (0.1322) 

0.1378 (0.1370) 

98.8 (109.4) 

125.2 (141.1) 

138.3 (151.4) 

0.1239 (0.1223) 
0.1252 (0.1238) 

0.1258 (0.1242) 

o(ECH) 

*121 

0.1320 (0.1320) 

129.9 (140.4) 

0.1254 (0.1237) 

Table 2 

Predictions of various approximants for the ratio R{s) (with n/ = 5) in the MS (c 2 = 1.475) and the ’t 
Hooft (c 2 =0) scheme (in brackets). y / Q 2 = 34 GeV; RScl chosen: 9 rsc1 = ~Q 2 ) a (~Q 2 ) =0.0452. 

Approx. 

d(~Q z ) 

C- rf[3i(-Q^s,,) pr - 

r pr - 

'3 

(TPS) 

0.048016 (0.047976) 

0.048016 (0.047976) 


mu 

[2/2]^ 2 

W' 

0.047996 (0.047974) 
0.047978 (0.047965) 

0.047931 (0.047939) 

-4.72 (-0.56) 0.047996 (0.047974) 

-8.48 (-2.24) 0.047981 (0.047967) 

-15.42 (-6.89) 0.047952 (0.047948) 

-93.91 (-89.75) 
-97.67 (-91.43) 

-104.61 (-96.08) 

J(ECH) 

m _ 

0.047959 (0.047959) 

-7.65 (-3.50) 0.047984 (0.047962) 

-96.84 (-92.69) 


4. Summary 

For a given SAi (TPS) of an observable S, we 
can construct an RScl-independent approximant 
which reproduces that TPS to the given order 
~a n+1 . It is global, i.e., it goes beyond the (poly¬ 
nomial) form of the TPS and could thus give us 
some clues to the nonperturbative effects - possi¬ 
bly in contrast to the local methods (ECH, PMS). 
It is, in principle, an improvement over another 
global approximant - the diagonal Pade approxi¬ 
mant (dPA), the latter being RScl-invariant only 
in the large- /3q limit. The question of how to elim¬ 
inate the second major source of unphysical de¬ 
pendence, the scheme (c 2 -)depenclence, remains 
open. One possibility would be to choose an “op¬ 
timal” c 2 (open problem). Another would be to 
extend the method so as to give us, in a global 
manner, approximants that are simultaneously 
RScl- and C 2 -independent (an open problem). 
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